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The effect of stellar-mass black holes on the structural 
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ABSTRACT 

We present the results of realistic iV-body modelling of massive star clusters in the 
Magellanic Clouds, aimed at investigating a dynamical origin for the radius-age trend 
observed in these systems. We find that stellar-mass black holes, formed in the su- 
pernova explosions of the most massive cluster stars, can constitute a dynamically 
important population. If a significant population is retained (here we assume complete 
retention), these objects rapidly form a dense core where interactions are common, 
resulting in the scattering of black holes into the cluster halo, and the ejection of black 
holes from the cluster. These two processes heat the stellar component, resulting in 
prolonged core expansion of a magnitude matching the observations. Significant core 
evolution is also observed in Magellanic Cloud clusters at early times. We find that 
this does not result from the action of black holes, but can be reproduced by the effects 
of mass-loss due to rapid stellar evolution in a primordially mass segregated cluster. 

Key words: stellar dynamics - globular clusters: general ~ methods: A^-body simu- 
lations - Magellanic Clouds. 



1 INTRODUCTION 
> 

. ^ Globular clusters are central to a wide variety of astrophysi- 
. cal research, ranging from star formation, stellar and binary 
$_( ' star evolution, and stellar dynamics, through to galaxy for- 
, mation and evolution, and cosmology. These objects there- 
fore constitute an integral part of our understanding of the 
Universe, and it is clearly vital that their internal evolu- 
tionary processes are well understood. The Galactic glob- 
ular clusters, while close, are exclusively ancient objects 
(r > 10^" yr). We can therefore accurately determine the 
end-points of their evolution, but must infer the complete 
long-term development which brought them to these ob- 
served states. To directly observe cluster evolution, we must 
switch our attention to the Magellanic Clouds (LMC/SMC), 
which both possess extensive systems of star clusters with 
masses comparable to the Galactic globulars, but crucially 
of all ages: 10^ < r < 10^° yr. These systems are of funda- 
mental importance because they are the nearest places we 
can observe snapshots of all phases of cluster development. 

bison. Freeman fc Lauej (igsi), discovered a striking 
relationship between core radius [tc) and age for LMC clus- 
ters - namely that the observed spread in increases dra- 
matically with increasing age. Here, Tc is the observational 
core radius, defined as the projected radius at which the 
surface b rightness has decreased to half its central value. 
Recently, iMackev fc Gilmord l|2003al lbh used Hubble Space 



Telescope (HST) WFPC2 imaging of 63 massive Magellanic 
Cloud clusters to more clearly demonstrate the radius-age 
trend in the LMC and show, for the first time, that a radius- 
age trend also exists in the SMC. An additional 46 objects 
have since been observed with HST/ACS (Program #9891) 
to improve sampling of the radius-age plane. Structural mea- 
surements for all 107 clusters may be seen in Fig[T] 

The observed radius-age relationship provides strong 
evidence that our understanding of globular cluster evolu- 
tion is incomplete, since standard quasi-equilibrium models 
do not predict large-scal e core expansion spann ing a full 
cluster life-time (see e.g.. lMevlan fc Heggielll997l ). Discern- 
ing the origin of the radius-age trend is therefore of consid- 
erable importance. A number of groups have investigated 
ossible explanation s - these include a size-of-sample bias 
Hunter et al.ll2003l). heating due to binary stars or tidal 
shocks ( Wilkinson et al.ll2003l ). and the formation of cores 
in primordially cusped clusters duo to the sinking of mas- 
sive stellar remnants (Mcrritt ct al. 2004). However, a model 
which fully accounts for the observed trend has yet to be 
elucidated. 

The radius-age trend is indistinguishable in the LMC 
and SMC, and the oldest LMC/SMC clusters have rc dis- 
tributions consistent with those of globular clusters in our 
Galaxy and in the Fornax and Sagittarius dSph galaxies 
(|Mackev fc Gilmord l2003 l. Since these galaxies have very 
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different tidal fields and possible external torques, this sug- 
gests that the radius-age relation is driven by internal clus- 
ter processe s, with any externa l or tidal effects second order 
(see also .Wilkinson et al.llioosi '). In this Letter, we report on 
the results of direct, realistic A'^-body simulations designed 
to investigate an internal dynamical origin for the radius- 
age trend. We follow the structural evolution of model clus- 
ters with varying degrees of primordial mass segregation 
(MSeg), possessing populations of stellar-mass black holes 
(BHs) formed in the supernova explosions of the most mas- 
sive cluster stars. We demonstrate that a cluster which re- 
tains its BHs undergoes dramatic core expansion for most 
of its lifetime, in contrast to a cluster with no BHs, which 
proceeds towards core collapse. We also show that primor- 
dial MSeg has an important effect on the early evolution of 
a cluster, when mass-loss due to stellar evolution is severe. 

2 NUMERICAL SETUP 

Direct A''-body modelling is a powerful tool for studying star 
cluster evolution because it incorporates all relevant physics 
with a minimum of simplifying assumptions. We have used 
the NBODY4 code in comb ination with a 32-ch ip GRAPE-6 
special-purpose computer (jMakino et ahlfioOSl ) to run simu- 
lations of Magel lanic Cloud clu sters. Full details of nbody4 
are provided bv I AarsethI l|2003l ). It uses a fourth-order Her- 
mite scheme and evaluation of the force and its first time 
derivative by the GRAPE-6 to integrate the equations of 
motion. Close encounters between stars, including stable bi- 
nary systems, are treated with two-body or chain regulariza- 
tion algorithms. Also incorporated are routines for modelling 
the s t ellar e volution of single and binary stars l|Hurlev et al.l 
I2OOOI . I2OO2I '). These include a metallicity dependence, and 
a mass-loss prescription such that evolving stars lose gas 
through winds and supernova explosions. 

We generate models with initial properties as close as 
possible to those observed for young Magellanic Cloud clus- 
ters. These objects possess radial surface brightne ss (SB) 
profiles best described by lElson. Fall fc FreemanI l| 19871 . 
hereafter EPF) models: ^i{r) = /io(l -I- r^/a^)""^/^, where 
/io is the central SB, a is the scale length, and 7 the 
power- law fall-off at large r. Typically, their core radii 
rc = a(2^/T - 1)^/^ ~ 0.2- 2.5 pc, and 7 - 2.0-3.5 
fe.g.. [Mackev fc Gilmorel l2003al ) . Their central densities po 
(M0pc~^) lie in the range 1.5 < logpo 2.5 (except for 
R136 which is much denser with logpo ~ 4.8), while their 
total masses Mtot (Mq) lie in the rang e 4.0 < log Mtot 
5.6 (McLaug hlin fc van der Marelir2005i ). We generate non- 
MSeg clusters by selecting stellar positions randomly from 
the density distribution of an EPF model with 7 = 3. Each 
star is assigned a velocity drawn from a Maxwellian distri- 
bution, where the velocity dispersion a is calculated using 
the Jeans equations assuming an isotropic velocity distri- 
bution. Expressions for a are given in a forthcoming paper 
(Mackey et al. 20 07a, in prep.). We select the initial mass 
function (IMF) of iKroupa (|200lh . with a stellar mass range 
0.1 — 100 M0. Choosing A'^ ~ 10^ particles results in cluster 
masses of log Mtot ~ 4.75. We adopt [Pe/H] = 0, similar to 
young LMC clusters. However, both Clouds exhibit strong 
age- metallicity relationships - this may have important im- 
plications for our results. 

Resolved observations of very young Magellanic Cloud 



clusters invariably reveal some degree of MSeg (e.g., 
Ide Griis et al.|[200a ). Detailed cl uster formation models sup- 
port such observations (e.g., iBonnell fc Bate! |2006| ): we 
would therefore like to include the effects of MSeg in our 
modelling. We have developed a method to generate clus- 
ters with primordial MSeg in a "self-consistent" fashion; 
again, full details will be provided by Mackey et al. (2007a, in 
prep.). Briefly, we take a non-MSeg cluster and use nbody4 
to evolve it without stellar evolution, so that the cluster be- 
gins to dynamically relax. The degree of primordial MSeg is 
controlled via the duration of this "pre-evolution" , Tms- The 
positions and velocities of the stars in the pre-evolved clus- 
ter are then used as the initial conditions (r = 0) for a full 
run with stellar evolution included. Stars slowly escape dur- 
ing pre-evolution, so our MSeg models are marginally less 
massive than non-MSeg models. We generate MSeg clusters 
with Tms = 450 Myr. These models have structural proper- 
ties (e.g., density profile, and radial mass-function variation) 
consistent with those observed for very young Magellanic 
Cloud clusters. 

LMC clusters are observed at galactocentric radii span- 
ning ~ — 14 kpc; our models move on circular orbits of ra- 
dius 6 kpc about a poin t-mass LMC with Mg = 9 X 10^ M©. 
IWilkinson et al.l (|2003l ) describe the implementation of an 
external tidal field within nbody4. Adopt ing a point-mass 
LMC is an over-simplification; however I Wilkinson et al.l 
(2003) showed that a weak external field does not result 
in strong core evolution - hence, here we are only interested 
in internal processes. Clusters are assumed to initially just 
fill their tidal radii. The initial tidal radius of a model clus- 
ter therefore set s the rat i o bet ween the length units used 
by NBODY4 (see lAarsethl |2003| '1 and physical length units 
(pc). This scaling controls the physical density of the clus- 
ter and hence the physical time-scale on which internal dy- 
namical processes occur. Our non-MSeg clusters have central 
density logpo = 2.31 and core radius rc ~ 1.90 pc, which 
matches typical young LMC and SMC objects. The primor- 
dially MSeg clusters have logpo = 4.58 and rc — 0.25 pc, 
which closely resembles the compact, massive LMC cluster 
R136. Given this correspondence, we are confident in our 
selection of an appropriate length-scale. 

We have modified nbody4 to control the production 
of BHs in supernova explosions. We can vary the minimum 
mass of a BH progenitor star, the masses of the BHs them- 
selves, and the natal velocity kicks they receive. This is im- 
plemented in a simple but serviceable manner. All stars ini- 
tially above 20 Mq produce BHs, with masses uniformly dis- 
tributed in the range 8 ^ rriBH ^ 12 Mq. This range is con- 
sistent with dynamical masses obtained from observations of 
X-ray binaries (e.g.. iCasares 2006). Each model cluster has 
the same random seed and so each begins with an identical 
stellar population: our adopted IMF and total A'' lead to the 
formation of 198 BHs in all clusters. In our models, natal BH 
kicks are either much larger than the cluster escape velocity 
Vosc (i.e., BH retention fraction /bh = 0) or zero (/bh = 1). 

To obtain structural measurements consistent with 
those for real clusters, we simulate observations of our A''- 
body models. That is, we mimic the reduction procedures 
from which the HST Vc measurements were derived. In those 
observations, the bright (saturation) and faint (background- 
limited) stellar detection levels are a weak function of cluster 
age, reflecting the requirement for longer exposure durations 
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Table 1. Details of A'^-body runs and initial conditions. Each cluster begins with A'^o stars with masses summing to Mtot, and 
initial central density po. Initial cluster structure is "observed" to obtain Tc and 7. Each model is evolved until Tmax- 



Name 




log A/tot 


log Po 




7 


Initial MSeg 


BH Retention 








(M©) 


(Mqpc-S) 


(pc) 




(Tms) 


(/bh) 


(Myr) 


Run 1 


100 881 


4.746 


2.31 


1.90 ±0.09 


2.96 ±0.17 


None 


0.0 


12 000 


Run 2 


100 881 


4.746 


2.31 


1.90 ±0.09 


2.96 ±0.17 


None 


1.0 


10 668 


Run 3 


95 315 


4.728 


4.58 


0.25 ±0.04 


2.33 ±0.10 


450 Myr 


0.0 


11274 


Run 4 


95 315 


4.728 


4.58 


0.25 ±0.04 


2.33 ±0.10 


450 Myr 


1.0 


10 000 



to image main sequence stars in older clusters. This is not 
responsible for the radius-age trend, but must be accounted 
for in our analysis. Further, the WFPC2 and ACS fields-of- 
view limit radial profiles to a maximum extent of ~ 70 — 100 
arcsec. To simulate these observations we first convert the 
luminosity and effective temperature of each A^-body star 
to magnitude and colour using the mod el atmospheres of 
iKurucj j 19921 ') and lBergeron et all (llQQSl ). Next, we impose 
appropriate bright and faint detection limits along with the 
field-of-view limits. We use the remaining s t ars to c onstruct 
a SB profile, following iMackev fc Gilmord l|2003ah . Stellar 
positions are projected onto a plane, and the SB calculated 
in circular annuli about the cluster centre. A varying annu- 
lus width is used to evenly sample both the cluster core and 
halo. Finally, we fit an EFF model to the resulting profile to 
derive and 7. To reduce noise we average the results for 
three orthogonal planar projections. 

3 SIMULATIONS AND RESULTS 

The parameter space of interest is spanned by non-MSeg 
clusters and those with significant primordial MSeg. For 
each of these, we consider evolution with BH retention frac- 
tions /bh ~ (large natal kicks) and 1 (no natal kicks). 
These four runs define the extremities of the parameter 
space, and hence are expected to cover the limits of cluster 
behaviour. Their properties are listed in Table [T] No spe- 
cial significance should be attached to Tmax - these simply 
represent the most convenient termination points for each 
simulation after r = 10 Gyr had been reached. 

3.1 A-body pair 1: No mass segregation 

The evolution of our non-MSeg runs is illustrated in Fig. [T] 
Run 1 constitutes the simplest case - no primordial MSeg 
and no retained BIfs. It behaves exactly as expected for a 
classical globular cluster. There is an early mass-loss phase 
(''' ^ 100 Myr) due to the evolution of the most massive 
cluster stars. During this phase, BHs are formed in super- 
nova explosions between 3.5-10 Myr; however, all receive 
large velocity kicks and escape. The severe early mass-loss 
is not reflected in the evolution of rc, presumably because it 
is evenly distributed throughout the cluster. Subsequently, 
the core radius slowly contracts as two-body relaxation pro- 
ceeds and mass segregation sets in. The median relaxation 
time at T = 10* yr is Uh ^ 2 Gyr. At Tmax ~ 12 Gyr « 6 trh 
the cluster has not yet entered the core-collapse phase. 

Now consider Run 2, which is identical to Run 1 ex- 
cept that /bh = 1. Once early stellar evolution is complete, 
the BHs are more massive than all other cluster members 
(of mean mass m, « 0.5 M©) and are hence subject to mass 
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Figure 1. Core radius evolution of Runs 1 and 2, which are 
initially identical, with no primordial MSeg. They have /bh = 
and 1, respectively. Run 1 evolves exactly as expected, gradually 
contracting as it moves towards core collapse. In contrast, the BHs 
in Run 2 induce dramatic core expansion after Ri 650 Myr. The 
plotted LMC/SMC data consists of all clusters from the WFPC2 
study o f lMackev fc Gilmorj l|2003al lb[) as well as preliminary ACS 
results from Mackey et al. (2007b, in prep.) (program #9891). 

stratification on a time-scale of ~ (m«/mBH) trh ~ 100 Myr. 
By 200 Myr, the mass density of the BHs within a radius of 
0.5 pc is already roughly equal to that of the stars; by 400 
Myr it is about three times larger. Soon after, the cen tral BH 
subsystem becomes unstable to further contraction l|Spitzeij 
Il987 , Eq. 3-55) and decouples from the stellar core in a run- 
away gravothermal collapse. At 490 Myr, the central density 
of the BH subsystem is ~ 80 times that of the stars. This 
is sufficient for the creation of stable BH binaries in three- 
body interactions - the first is formed at ~ 510 Myr, and 
by 800 Myr there are four. Until this phase, the evolution of 
Run 2 is observationally identical to that of Run 1. Neither 
BH retention, nor the subsequent formation of a central BH 
subsystem leads to differential evo l ution of Vc- This contrasts 
with the results of iMerritt et al.l l|2004h who found signifi- 
cant early expansion in their models due to the sinking of 
BHs. We attribute this difference to the much higher degree 
of central mass concentration in their initially cusped clus- 
ters, which thereby respond more strongly and more rapidly 
to the perturbations induced by sinking remnants. These 
authors also noted the possibility of further cluster expan- 
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sion due to subsequent evolution of the BH subsystem. We 
indeed observe expansion due to such processes (see below) . 

Once formed, binary BHs undergo superelastic colli- 
sions with other BHs in the core. The binaries become 
"harder" , and the released binding energy is carried off by 
the interacting BHs. This leads to BHs being scattered out- 
side rc, often into the cluster halo, as well as to BHs being 
ejected from the cluster (we retain this terminology hence- 
forth). Eventually a BH binary is sufficiently hard that the 
recoil velocity imparted to it during a collision is larger than 
the cluster escape velocity, and the binary is ejected. A BH 
scattered outside the cluster core gradually sinks back into 
the centre via dynamical friction, thus transferring its newly- 
gained energy to the stellar component of the cluster. Most 
is deposited within where the stellar density is greatest. 
The ejection of BHs also transfers energy to the cluster, since 
a mass m escaping from a cluster potential well of depth |<l?| 
does work m|$| on the cluster. This mechanism is partic- 
ularly effective in heating the stellar core, since BHs are 
ejected from the very centre of the cluster, and the energy 
contributed to each part of the cluster is proportional to the 
contribution which that part makes to the central potential. 

Together, these two processes (scattering and ejection) 
result in significant core expansion, starting between r « 
600 — 700 Myr. Expansion continues for the remainder of the 
simulation, which terminates at Tmax ~ 10.6 Gyr. The size of 
rc is roughly proportional to log r, consistent with the upper 
envelope of the observed cluster distribution. However, in 
this model the expansion begins too late for the evolution 
to trace the upper envelope exactly; rather, it runs parallel. 

The number of stable BH binaries in the system peaks 
at 5 at T « 890 Myr. After this point, there are — 5 
BH binaries at any given time. Single and binary BHs are 
continually ejected; however, empirically, both escape rates 
depend logarithmically on r - i.e., dNc/dr oc logr. This 
arises due to the decreasing density of the central BH sub- 
system - the number of BHs is falling because of ejections; 
these ejections also heat the BH core. The BH-BH encounter 
rate therefore decreases with time. Hence, the BH binary 
hardening rate decreases, as do the BH ejection rates. Fur- 
thermore, the stellar core is also less efSciently heated with 
time - this is reflected in the roughly logarithmic depen- 
dence of Tc on T. By Tmax ~ 10.6 Gyr, 96 single BHs, 15 
binary BHs and one triple BH have escaped; 65 single BHs 
and 2 binary BHs remai n in the cluster. This is at odds 
with early studies fe.g.. iKulkarni. Hut fc McMiUanl 1 19931 : 
ISigurdsson fc Hernauistll 19931 ) which predicted depletion of 
BH populations on timescales much less than cluster life- 
times. The decreasing BH encounter rate seen in our models 
prolongs the life of the BH subsystem for much longer than 
previously appreciated. 

The mean mass of stellar escapers is identical in both 
Runs 1 and 2, at 0.339 Mq. This is less than m, at all times. 
The distributions of velocities with which stars escape are 
also indistinguishable. These results imply that both mod- 
els lose stars solely due to relaxation processes. There is 
only a tiny group of ~ 30 high velocity escapers in Run 
2, indicating that stars interact closely with BH binaries 
only very rarely. Heating of the stellar component via close 
interactions is negligible - the hardening of BH binaries is 
driven solely through interactions with other BHs. At r = 10 
Gyr, the masses of Runs 1 and 2, respectively, are 0.36 Mtot 



and 0.29 Mtot, reflecting the fact that Run 2 is more loosely 
bound than Run 1 for the majority of its evolution. 

3.2 TV-body pair 2: Strong mass segregation 

Runs 3 and 4 are primordially MSeg versions of Runs 1 and 
2, respectively. Early mass-loss due to stellar evolution is 
highly centrally concentrated - hence the amount of heat- 
ing per unit mass lost is maximised, leading to dramatic 
early core expansion (Fig. [2]). Run 3 traces the observed up- 
per envelope of clusters until several hundred Myr. Run 4 
retains its BHs and hence loses less mass than Run 3 - this 
is reflected in its smaller rc- After the early mass- loss phase 
is complete, core expansion stalls in both runs. Two-body 
relaxation gradually takes over in Run 3, leading to a slow 
contraction in Vc- At r = 1 Gyr, trh ~ 4 Gyr; hence this 
cluster is not near core collapse by Tmax ~ 11.2 Gyr. At 
r = 10 Gyr, the remaining mass in Run 3 is 0.30 Mtot. 

In Run 4, the BH population evolves similarly to that 
in Run 2. One might naively expect the earlier development 
of a compact BH subsystem in Run 4, because the BHs 
are already located in the core due to the primordial MSeg. 
However, the centrally concentrated mass-loss hampers the 
accumulation of a dense BH core, and the first binary BH 
does not form until 570 Myr, a similar time to the non-MSeg 
model. The BH subsystem evolves more slowly than that in 
Run 2 - by Tmax ~ 10 Gyr, there are still 95 single BHs 
and 2 binary BHs remaining in the cluster. As in Run 2, 
the evolution of the BH subsystem leads to expansion of Vc- 
This begins at r ~ 800 Myr and continues until Tmax- As 
previously, Vc behaves roughly as log r during this phase. 
By Tmax, Run 4 has rc ~ 11 pc, comparable to that ob- 
served for the most extended old Magellanic Cloud clusters 
(e.g.. Reticulum). However, it is only weakly bound, retain- 
ing ~ 0.13 Mtot- This mass loss is not the driver for the core 
expansion, so a more massive cluster could have compara- 
ble expansion while retaining more of its total mass. Indeed, 
extended old LMC clusters typically have masses ~ 
(Mackcv fc Gilmore|[2003al ). which may easily be ~ 0.1 Mtot- 

4 DISCUSSION 

Our four simulations cover the observed cluster distribution 
in radius-age space, thereby defining a dynamical origin for 
the radius-age trend. At ages less than a few hundred Myr, 
cluster cores expand due to centrally concentrated mass-loss 
from stellar evolution. At later times, expansion is induced 
via heating due to a BH population. Although early mass- 
loss may result in significant core expansion, a cluster can- 
not reach the upper right corner of the radius-age space by 
this means alone - the mass-loss phase is too short, and the 
maximum allowed expansion rate during this phase is de- 
fined by the observed upper envelope of clusters. Only with 
prolonged expansion due to BHs can ~ 10 Gyr old clusters 
with rc > 6 pc be explained in this model. Although we have 
assumed /bh = 1, full retention is not necessary for cluster 
expansion. BH kicks of order 10 < Wkick ^ 20 kms~^ would 
result in /bh ~ 0.5 in our models; we expect rc evolution 
in such systems to be intermediate between that of Runs 1 
and 2, or Runs 3 and 4. We will address this issue further 
in an upcoming paper (Mackey et al. 2007a, in prep.). 

Galactic globular clusters, with TV ~ 10*', are an order of 
magnitude more massive than our present models. However, 



Black holes in massive star clusters 5 




logio (t / Myr) 



Figure 2. Core radius evolution of Runs 3 and 4, which are 
initially identical, with significant primordial MSeg. They have 
/bh = E'lid 1, respectively. Both expand dramatically at early 
times due to mass-loss from stellar evolution. Subsequently, Run 
3 begins to contract as two-body relaxation proceeds. In contrast, 
Run 4 continues expanding due to its BH population. 

we expect the evolution described above to scale to such 
objects - reflected in the fact they possess an rc distribu- 
tion consistent with that observed for the oldest Magellanic 
Cloud clusters (|Mackev fc Gilmord liooj ). This is because 
the mass fraction of BHs formed in a cluster is dependent 
only on the IMF and minimum progenitor mass, neither of 
which should change with Mtot , while a larger A/tot implies a 
larger /bh since it is easier to retain newly-formed BHs. The 
densities in our models are consistent with those observed 
for globular clusters; hence we expect the same processes to 
operate on similar time-scales, although BHs are likely to be 
more difficult to eject in more massive clusters - increasing 
the potential of each BH to heat the cluster via additional 
scattering-sinking cycles. Core expansion due to mass-loss or 
BH heating has strong implications for the observed proper- 
ties of Galactic globular clusters (e.g., the fraction which are 
core-collapsed) as well as their survivability. Extended clus- 
ters are significantly more susceptible to tidal disruption, so 
it is important to account for expansion effects in studies 
of the evolution of the globular cluster mass function, for 
example. Core expansion due to BHs may also offer a vi- 
able explanation for the origin of the luminous, unusually 
extended globular clusters f ound in M31, whic h are > 10 
Gyr old metal-poor objects l|Mackev et al.|[2006l) . 

Our model requires variations in BH population size 
between otherwise similar clusters. There are a variety of 
possibilities in this regard. First, the number of BH- forming 
stars in a cluster is small, so there will be sampling- noise 
variations between clusters. Further, any dispersion in stel- 
lar rotation may introduce mass-loss variations and fur- 
ther dispersion in BH numbers. Natal BH kicks are poorly 
constrained at present - typical estimates lie in the range 
J; Wkick ^ 200 kms~^ , with kicks of a few tens of kms~^ 
possibly favoured fe.g.. IWillems et aL I l2005l . and references 
therein). Stellar binarity may therefore play a significant 



role in retaining cluster BHs, as will the initial cluster mass 
and degree of primordial MSeg, especially if Hkick ~ «csc. 
Metallicity may also be a key factor, as theory suggests that 
BH production is more frequent, and ttibh i s greater for 
meta l poor stars than for metal rich stars (e.g.. IZhang et aU 
l2007t ). In this respect, the age-metallicity relationships of 
the Magellanic Clouds (where [Fe/H] decreases for clusters 
of increasing age) may play a central role in shaping the 
radius-age trend. Similarly, the spread in [Fe/H] for Galac- 
tic globulars may have been important in determining the 
structural properties of these objects. Our results imply that 
clusters possessing significant BH populations are, for most 
of their lives, low-density objects in which the timescale for 
close encounters between stars and BHs is very long. It is 
therefore unsurprising that no BH X-ray binaries are s een 
in the ~ 150 Galactic globulars l|Verbunt fc Lewinllioosl '). 
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